(3 + l)-spectrum of neutrino masses: A chance for 

LSND? 



O. L. G. Peres^'^ Q and A. Yu. Smirnov^'^ 

(1) The Abdus Salam International Centre for Theoretical Physics, I-S^lOO Trieste, Italy 

(2) Instituto de Fisica Gleh Wataghin, Universidade Estadual de Campinas, UNICAMP 
§ ■ 13083-970 Campinas SP, Brazil 



(N 

(N 
(N 



(3) Institute for Nuclear Research of Russian Academy of Sciences, Moscow 117312, 

Russia 



Abstract 



> 

' If active to active neutrino transitions are dominant modes of the atmospheric 

(f^ — > Ur) and the solar neutrino oscillations (fg ^fi/^r), as is indicated by recent 
I data, the favoured scheme which accommodates the LSND result - the so called (2+2)- 

(—] ' scheme - should be discarded. We introduce the parameters t^^**" and r/l"" which 

I ' quantify an involvement of the sterile component in the solar and atmospheric neutrino 

^ • oscillations. The (2 + 2)-scheme predicts r/"*™" + r/l"" = 1 and the experimental proof of 

deviation from this equality will discriminate the scheme. In this connection the (3+1)- 
^ ■ scheme is revisited in which the fourth (predominantly sterile) neutrino is isolated from 

H ; 

5^ I a block of three flavour neutrinos by the mass gap AmlsND ~ (0.4-10) eV^. We find 

that in the (3 + l)-scheme the LSND result can be reconciled with existing bounds on 
fg- and - disappearance at 95—99% C.L.. The generic prediction of the scheme is the 
Ve- and - disappearance probabilities at the level of present experimental bounds. 
The possibility to strengthen the bound on v^- disappearance in the KEK - front 
detector experiment is studied. We consider phenomenology of the (3 + l)-scheme, 
in particular, its implications for the atmospheric neutrinos, neutrinoless double beta 
decay searches, supernova neutrinos and primordial nucleosynthesis. 
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1 Introduction 



It is widely accepted that simultaneous explanation of the solar 0,11, atmospheric 
and LSND results in terms of neutrino conversion requires an existence of the sterile 
state Less appreciated and discussed fact is that explanation of the LSND result in the 
favored 4z/-schemes implies that the solar or atmospheric neutrinos (or both) are converted 
to sterile state. Indeed, it is claimed that the only 4i/-schemes which can explain the 
LSND result and accommodate oscillation solutions of the solar and atmospheric neutrino 
problems are the schemes of the (2 + 2) type. According to the (2 + 2) - scheme, four 
neutrino mass eigenstates form two nearly degenerate pairs separated by the gap Am^ = 
Am|5jvD ~ 1 eV'^. Splittings between masses within pairs are much smaller: they are 
determined by Am^^^ ^ 3- 10~^ eV^ in the atmospheric neutrino pair and Am| < 10"^ eV^ 
in the solar neutrino pair. 

The key feature of the (2 + 2)-scheme which allows one to get a large enough probability 
of the z/^ — > Ue transition and to satisfy existing bounds on mixing parameters is that z/g and 
z/^ are present in different pairs: Ue is mainly distributed in the pair with Ami, - splitting, 
whereas z/^ is mainly distributed in the pair with Am^^^ - splitting. Indeed, the depth of 
the z/^ ^ z/e oscillations driven by Am|^^^ equals 

sin^ 2eLSND = 4 ( ^ U^,U,,y = 4(5] U^.U,,)' ~ e', (1) 

i=3,4 j=l,2 

where summation is over the mass eigenstates in the heavy (or light) degenerate pair. If we 
introduce small parameter e which describes admixture of z/g (or u^) in the pair where it is 
not a dominant component, then taking into account that in a given pair either U^j ~ 0(1) 
or Uej ~ 0(1)) "we get that the effective mixing is of the order e^. 

The sterile neutrino can be distributed in two pairs in different ways, and there are two 
extreme versions of the (2 + 2)-scheme. 1) The sterile neutrino is mainly in the solar pair, 
so that the z/g Vs conversion is responsible for the solution of the solar neutrino problem, 
whereas atmospheric neutrino problem is solved by z/^ ^ oscillations. 2) The sterile 
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neutrino is in the atmospheric pair, then z/^ oscillates into z/^ and the solar neutrino problem 
is solved by the ^ conversion. The intermediate situations are possible when both the 
solar and atmospheric neutrinos partly oscillate to the sterile component. The exceptional 
case is when is distributed equally in the solar and atmospheric neutrino pairs. In all 
other cases the sterile neutrino contributes either to solar or to atmospheric neutrinos more 
than one half. 

In the alternative (3 + l)-scheme ( see Fig. |1]), three mass eigenstates form "flavour 
block" which has predominantly active flavour with small admixture of the sterile state and 
small mass splitting. The fourth mass eigenstate is separated from the flavour block by 
the LSND mass gap Am|^^jy^ (0.4 — 10) eV^; it consists mainly of the sterile neutrino 
with small admixtures of active neutrinos: \Usi\^ ~ 1, |f/Q4p << 1, a = e,/i and r. (In 
general, there is mixing of Vg and Vr in ^^4 and in the flavor block [jl^].) In the (3 + l)-scheme 
the atmospheric and solar neutrino data are explained by the oscillations and 

—>■ i^fj./ conversion respectively. 

It is claimed that the (3 + l)-mass scheme can not reproduce large enough probability 
of the z/^ Ve oscillations to explain the LSND result Indeed, the depth of z/^ 

oscillations driven by the Am^^^^, equals 

sin2 2^^,.5iVD = 4f/^24^7i~e^ (2) 

where and [/e45 are the admixtures of the v^. and in the 4th mass eigenstate. These 
admixtures are small (~ e) being restricted by the accelerator and reactor experiments. 
As a consequence, the depth of oscillations is of the fourth order in the small parameter. 
Notice, however, that early analysis of data (with small statistics) showed that the bounds 
from accelerator and reactors experiments and the positive LSND result can be reconciled 
in the context of one level dominance scheme at 99% CL. ||11|| . 

There are three recent results which may change eventually the situation in favor of the 
(3 + l)-scheme: 
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1. The data on atmospheric neutrinos (specifically on zenith angle distribution of the 
upward-going muons, on the partially contained multi-ring events, with > 5 GeV, 
and on the enriched neutral current event sample) disfavour oscillations to the ster- 



ile state [|T2|. The ^ Ug oscillations can be accepted by the data at 3a level. 



Oscillations of active neutrinos, i/^ — > z/,-, give better fit of the experimental results. 

2. Conversion of solar z/g neutrinos to active neutrinos gives better global fit of experi- 
mental results than conversion to sterile neutrinos 0. 

If active neutrino channels dominate both in the solar and in atmospheric neutrino 
oscillations, the (2 + 2) -mass scheme should be discarded and the oscillation in- 
terpretation of the LSND result becomes problematic. In this connection, we have 



reconsidered the possibility to explain LSND results in the (3 + 1) -mass scheme [O 



3. Latest analysis of the LSND data 0] shows in further shift of the allowed region 
of oscillation parameters to smaller values of mixing angle. This leads to a better 
agreement between the bounds obtained in CDHS CCFR [|1^] and Bugey |T6 



experiments and the LSND result in the context of (3 + l)-mass scheme |10, 17, IS . 



In this paper we elaborate on the phenomenology of the (3 + l)-scheme which can 
accommodate the oscillation interpretation of the LSND result. In sect. 2 we introduce 
parameters r^^"" and 77^*"^ which quantify participation of the sterile neutrino in the solar 
and atmospheric neutrino oscillations. We show how these parameters can be used to 
disentangle the (2 + 2)- and (3 + l)-schemes. In sect. 3 we analyze the u^j, z/g oscillation 
probability in the (3 + l)-scheme. We estimate the confidence level for the bound on sin^ 29efj, 
from z/^-and z/g-disappearance experiments. In sect. 4, we study a possibility to improve the 
bounds on z/^-disappearance, and consequently, on sin^ 26^6^ with the KEK - front detectors 
experiments. In sect. 5, we describe phenomenological implications of the (3 + l)-scheme for 
atmospheric neutrinos, neutrinoless double beta decay searches, primordial nucleosynthesis 
and supernova neutrinos. In sect. 6 we summarize our results and discuss perspectives to 
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identify (3 + l)-mass scheme. In Appendix, we evaluate how strongfy the LSND oscillation 
probability can be enhanced in the schemes with more than one sterile neutrino. 



2 The (2 +2)-scheme versus (3 + l)-scheme 

As we have pointed out in the introduction, the present data disfavor pure oscil- 
lations of the atmospheric neutrinos and moreover pure z/g Ug conversion is not the best 
solution of the solar neutrino problem. Let us consider a general case when Us is involved 
both in the solar and in the atmospheric neutrino oscillations. 



2.1 Excluding the (2 +2)-scheme ? 

For definiteness we will consider the (2 + 2)-scheme with z/i and i>2 being the "solar" pair of 
the eigenstates with AttIq = m\ — m\ mass splitting and and 1/4 - the atmospheric pairs 
(with Am^j^ = m\ — m| mass splitting (mi < m2 < < 1714). To a good approximation 
(as far as solar and atmospheric neutrinos are concerned) we can neglect small admixtures 
of Uf, and implied by the LSND result. So, the i/g-flavor is present in the solar pair: 

|t/eir + |t/e2|'~l, (3) 

and the z/^ flavor is in the atmospheric pair: 

\U^3\' + \U,a\'-1. (4) 

Under conditions (^) and (^ the effect of the sterile neutrino both in solar and in 
atmospheric neutrinos is described by a unique parameter [^. In what follows, we will 
introduce this parameter in a more transparent way. Notice that since both solar and 
atmospheric neutrino oscillations are reduced to 2z/ cases, there is no CP-violation effect 
and the absolute values of mixings are relevant only. 
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In general, Vr and Ug are mixed both in the solar and atmospheric neutrino pairs. Using 
the unitarity of the mixing matrix and the equality it is easy to show that in the solar 
pair (i^i, V2) the electron neutrino mixes with the combination: 

z> = cos a Vs + sin a Ur, (5) 

(where a is the arbitrary mixing angle) so that 

z/i = cos ^0 z/g — sin ^0 i>, 1^2 = cos 9q u + sin 9q . (6) 

Here ^0 is the angle which appears in the 2i/-analysis of the solar neutrino data. Then as 
follows from the unitarity and the condition (^), the combination of and orthogonal 
to that in Eq. d): 

V = cos a Ur — sin a Us, (7) 
mixes with z/^ in the atmospheric neutrino pair (1/3, z/4): 

Z/3 = cos 9atm^^, - sin OatmJ^', = COS 9atmJ^' + slu 9atmT^^„ (8) 

where 9atm, ~ 45*^ is the mixing angle responsible for oscillations of the atmospheric neutrinos. 
Using Eqs. (||, |^) we find 

rjl-n^P^^^2^P^^^2^^^^2^^ (9) 

that is, cos^ a or t]!"" give the total presence of the sterile neutrino in the solar pair and 
describe the effects of the sterile component in the solar neutrinos. Similarly, from the Eq. 
(0, 1) we get 

7^f'" = |t/,3|2 + |f/,4|2 = sin2a, (10) 
so, sin^ a or 77°*"* give a total presence of the sterile component in the atmospheric pair. 
Clearly, in the (2 + 2)-scheme: 

r/, = r^f " + T^f = 1. (11) 
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The schemes with pure z/g Vs conversion of the solar neutrinos and Vr oscillations 

of the atmospheric neutrinos correspond to cos^ a = 1. The opposite situation: the solar 
z/g — Vt and the atmospheric z/„ — Vs transitions corresponds to cos^ a = 0. 



Notice that in notations of Ref. [19 



cos a = C23 ■ C24. 

The equality (pUf ) is a generic property of the (2 + 2)-scheme and it can be checked 
experimentally. One should measure or restrict r^f"" from the solar neutrino data, and 
independently, r/^*™ from the atmospheric neutrino data. According to Eqs. and (^, 
y/r]f^ is the admixture of the Vg in the state 

y = y^rff^ Vs + y^l - Vr (12) 
to which the solar convert. Thus, rjf^^ should be determined from the fit of the solar 

.atm 



neutrino data in terms of the 2z/ mixing of and z/. Similarly (see Eqs. (|^), (|T0])), ^jrj] 
is the admixture of the z/, in the state 



^' = \/vf"' ^r- yl- vf"" (13) 

to which the atmospheric neutrinos oscillate. Thus, 77^*™ can be determined from the 2z/- 
analysis of the atmospheric neutrino data in terms of z/^ — u' mixing. 

If it will be proven that 

(or larger than 1), the (2 + 2)-scheme should be discarded. Let us summarize the present 
situation. The global analysis of the solar neutrino data P3| shows that for large mixing 
angle solutions (LMA): 

^sun ^ Q 44 (gg^^ Qj^^^ ^ ^sun ^ q (99% C.L.) ; (15) 

for LOW solution: 

r]f " < 0.30 (90% C.L.) , r/f " < 0.77 (99% C.L.) ; (16) 
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and for small mixing angle solution (SMA): 

r^f " < 0.90 (90%C.L.) , (17) 
and no bound appears at the 99%C.L.. The bound from the SK results on atmospheric 



neutrino data is ||35[| : 



r/f"^ < 0.67 (90%C.L.) . (18) 

One can get similar bound from the fit of the zenith angle distribution of events detected 
by MACRO m-. r/f" < 0.7 at 90 % C.L.. 



Thus, taking the LMA solution, we get 

r]s = vT" + < 1-11 (90%C.L.), (19) 

and for the SMA solution: r]s < 1.57 and for the LOW solution: rj^ < 0.97. That is, at 
the moment the (2 + 2)-scheme is well acceptable. However the forthcoming solar and 
atmospheric neutrino experiments can significantly strengthen this bound. 

Let us consider dependences of various observable on r^f"" and r/"*"^. 

The solar neutrinos undergo i) conversion. Difference of the Ug and z/ potentials in 
matter equals 

V = V2GFln,-VrY) ' ^20) 
where Gi;' is the Fermi coupling constant, rie and are the concentrations of electrons and 
neutrons correspondingly. With the increase of r^f™ the potential decreases, this leads to 
a shift of the adiabatic edge of the suppression pit to smaller Am?' / and to modification 
(weakening) of the Earth matter effect . 



One expects an intermediate situation between pure active and pure sterile cases. 

Presence of the sterile component in the solar neutrino flux modifies also interactions 
of neutrinos in detectors. The reduced rate of the neutral current events [NC] defined as 
the ratio of events with and without oscillations, Nnc{osc)/Nnc{SSM), decreases with 
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increase of t^™"-; 

here P is the effective ( averaged) survival probabihty. Two remarks are in order. 

1) . The probabihty P = P{rjf^'^) should be calculated with the effective potential (^OD , 
although its dependence on ry^"" is rather weak. 

2) . The probability should be calculated for each specific variable separately (taking 
into account cross sections, energy thresholds etc.) so that P in Eq. (pT]) may differ from 
P which will appear in the following formulas. 

The double ratio - the ratio of the reduced rates of the neutral current events [NC] and 
the charged current events ([CC] = Ncc{.osc)/Ncc{SSM)) changes as: 

[NC] _ _ 1 - r^r .221 

[CC] - + p ■ 

With increase of r/*"" the double ratio decreases from 1/P to 1 . The rate [NC] and the ratio 
[NC]/[CC] will be measured in the SNO experiment p? ]. 

The reduced charged current event rate [CC] (= Ncc{osc)/Ncc{SSM)) can be written 
in terms of the suppression factor for the z/ — e event rate, R^e = N^e{osc) / Ni^e{S S M) , as 

[CC] = ^ . (23) 

^ ^ l-r(l-r7f«")(l-l/P) ^ ^ 

Here r = cr(z/^e)/cr(z/ee) is the ratio of the cross-sections of the muon and the electron 
neutrino scattering on the electron. The rate R^e is measured with high precision at Super- 
Kamiokande and it will be also measured at SNO. According to (p3D , with increase of 
involvement of the sterile neutrino the rate increases from Rue/ (1~^(1~1/-P)) for the pure 
active case to R^e for the pure sterile neutrino case. 

The presence of the sterile component in the atmospheric neutrinos can be established 
by studies of the neutral current interaction rates. The suppression factor for the vr'^-event 



rate in the pure sterile case is expected to be about C, = 0.7 — 0.8 [38]. In general case the 
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ratio of vr'' to e-like event rates is suppressed as 

1 - + er/r ) , (24) 



where A^°(7r°) and are the numbers of events with and without oscillations. 

Appearance of the Vr in oscillations of is suppressed by the factor cos^ a = 1 — rjf"^. 
This can be tested in the long base-line experiments. One can compare the mixing parameter 
sin^ 29^fj_ extracted from studies of the z/^ - disappearance and the parameter sin^ 29^^ found 
from the Ur - appearance experiments. In the (2 + 2)-scheme one expects 

sin^ 29,, 



sin^ 29 



1 - 7]^"". (25) 



Also the zenith angle distribution of the of multi-ring events from the so called NC 
enriched sample [0] is sensitive to 1]° 



The Earth matter effect on atmospheric neutrino oscillations depends on r^"*™. The 
matter potential for z/^ — v' system equals: 

V = V2GFYVf""- (26) 
The potential increases with 77°*"^, and one expects an intermediate situation between active 



and sterile neutrino cases [p4| ,p5|. For high energy upward going muons the suppression 
of the oscillation effect increases with r^^*™-. The effect is more complicated for neutrinos 
crossing the core of the earth where the parametric enhancement of oscillations may take 
place. 



2.2 The (3 + l)-scheme and mixing of sterile neutrino 

In the (3 + l)-scheme with 4th (isolated) mass eigenstate being predominantly sterile one, 
total involvement of the sterile neutrino in the solar and atmospheric neutrino oscillations 
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can be very small. If the Vr component in the z/4 is absent, we have from unitarity condition 

3 

^ Ul = |f/e4|' + If/^r < (3 - 10) ■ 10-2 (27) 

i=l 

where the inequality corresponds to the CDHS and Bugey bounds. 



The situation is different if Vr mixes with in the 4th eigenstate |T^. Now 1/4 consists 
mainly of the combination 

z/ = cos [3 Us + sin [3 (28) 

with small admixtures of z/g and z/^ implied by the LSND experiment (z^4 ~ v). The mixing 
in the flavor block can be obtained from the one in the original (3 + l)-scheme (without 
Vs — mixing) by substituting z/^ — > z^', where 

v' = cos [3 Vr — sin [3 z/^. (29) 



is the orthogonal combination to that in Eq. (p8|). 

An immediate consequence of the Vr — Vs mixing is an appearance of Vr <-> Vs oscillations 
driven by Am^^^^, [0. The depth of these oscillations equals sin^ 2/5. 

Due to the presence of z/^ in z/4 described by f/^4 the scheme leads also to z/^ ^ Vr 
oscillations with Am^ = Am^^^j^ and the depth 

sin^ 20 ^ 4|?7^4|2 sin^ (3 = sin^ 29^^ ■ sin^ p , (30) 



where sin 26^^^ describes the disappearance of z/^. Eq. ( pO| ) allows to get the upper bound on 

sin^ p. Indeed, using sin^ 26 ^r{C / N) - the upper bound on sin^ 26^r from the CHORUS |^| 

.2 . 



and NOMAD |^ experiments and taking sin 26'^^ at the upper edge allowed by the CDHS 
experiment (as is implied by the LSND result) we get: 

From this formula we find that for Am^ < 2 eV^ no bound appears; for Am^ =4, 6, 8 
eV^ we get siv? (3 < 0.32, 0.13, 0.06 correspondingly. 
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In the presence of Vr — mixing the solar will convert into the combination 



z/^ cos Oatm + sin 6 atm COS /? - z/^ sin 6 atm sin /?. (32) 

If SMA is the solution for solar neutrino problem then the above state is roughly 2/2) . 

That is, the contribution of the sterile component is determined by the product sin Oatm ■ 
sin (3 and therefore 

r7r" = sin2^,i^sin2/3 (33) 
(instead of cos^a in the (2 + 2)-scheme). 

Due to Vj. — Ug mixing the atmospheric z/^ will oscillate with Am^ = Am^j^ into the 



state u' (see Eq. (^91) ). According to (|29|), the admixture of the sterile component is sin/3, 
so that 

r/f = sin^ (3. (34) 



From (|33|) and (^) we find that the total effect of the sterile components in the solar 



and atmospheric neutrino oscillations can be described by 

Vs = + VT = sin' /3 ■ (1 + sin^ Oatm)- (35) 

The parameter ris can be larger than 1 since in the (3 + l)-scheme we make the double 
counting of the Ug contributions from the two lightest states (they contribute both to t]!"" 
and to 7]^^). In contrast with the (2 + 2)-scheme, now 77^ is not fixed and it can change 
from (1 + sin^ Oatm) to practically zero when sin^ /3 ^ 0. 

According to Eqs. (131) and 



^-"-sin^^^i^r^f-, (36) 

so that the inequality 77^"" < 77^*™ is the generic property of the (3 + l)-scheme under 
consideration. Taking sin^ Oatm < 0.67 (which corresponds to the lower bound sin^ 26'ajm > 
0.88 from the Super-Kamiokande data) we get from ( ^3]) and (|18]): r/f"" < 0.45 independently 
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on solution of the solar neutrino problem. This bound is stronger than the immediate bound 
from the solar neutrino data. 



3 LSND oscillation probability in the (3 + l)-scheme 

Let us consider predictions for the — > oscillations relevant for the LSND result. 
The short range oscillation results are determined by one Am'j^g^j^ only (flavour block is 
"frozen") and the oscillation pattern is reduced to the 2z/ neutrino oscillations. So, one can 
use immediately the results of the analysis of the LSND data in terms of the 2z/-oscillations. 
The depth of the oscillations or the effective mixing parameter is given in Eq. (H). As in 
the one level dominance scheme [^], it is determined by the admixtures of and z/^ in the 



isolated state. This prediction does not depend on mixing in the flavour block, in particular, 
on the solution of the solar neutrino problem. 

Mixing matrix elements which enter the expression for sin^ 2^e^ (^) are restricted by 
short baseline experiments. In the range of Am^ relevant for the LSND result the best 
bound on U^^ is given by the Bugey reactor experiment |16[. The i/g-survival probability in 
this experiment is determined by 

sin^ 26,, = 4 U^^^il - U^J ^ 4 (37) 

The best bound on U"^^, for 0.3 eV^ < Am^ < 8 eV^, follows from the i/^-disappearance 



searches in the CDHS experiment WM. The relevant mixing parameter equals 



sin^ 29,, = 4 Ul,{l - Ul,) ^ 4 f/^,. (38) 

From Eqs. (0), (|37D and (^) we get unique relation between the depths of oscillations in 
the CDHS, Bugey and LSND experiments: 

sin^ 26,, ^ ^ sin^ 26,, ■ sin^ 26,,. (39) 
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This relation holds in the one level dominance scheme when only one Am^ contributes to 
oscillations in all three experiments. The relation is modified in more complicated situations, 
e.g., in schemes with more than one sterile neutrino. 



The Bugey [jT6[ and CDHS ||T^ experiments have published the 90% C.L. upper bounds 
on sin^ 26'ee and sin^26'^^ as functions of the Am^. These bounds transfer immediately to 
the bounds on U"^^ and U^^: 

<U'Bu,eyi^rn^), (40) 

<U^cDHsi^^')- 

Using Eq. (^) and Eq. (|38|) one gets the upper bound on sin^ 26^6^ 11)0: 

sin^ 26,, = 4 Uf,,U',, < 4 t/|„,,,(Am^) ■ U'aDHsi^rn'). (41) 



The bound ( |1T| ) shown in Fig. |^ excludes most of the region of parameters indicated by 
the LSND experiment. On this basis, it was concluded that (3 + 1) -mass scheme cannot 
reproduce the LSND result [|,|^. The question is, however, which confidence level should 
be ascribed to the bound (^)? Or, at which confidence level the LSND result is excluded 
by the combined Bugey and CDHS bounds? 

Let us introduce the probabilities, Pq, (f/^4, Am^), a = e, /x, that experimental data 
from Bugey or CDHS correspond to a given value of U^^ for fixed Am^. The probabilities 
Pa{U^4, Am^) should be found from the fit of experimental data. In such a fit ?7^4 can 
take both positive and negative values but with {U^^l < 1. Therefore the normalization 
conditions for Pa should be written as J^^dxPa{x, Am'^) = 1. The 90% C.L. bounds on 

are determined by the condition 

/ dxPa{x,Am^) = 0.9 . (42) 
J —00 

The probability that the product of U^^ ■ U^^ is larger than certain value p is given by 
the integral 

poc poo 

V{p)= dy dxPeix,Am^)P,iy,Am^). (43) 

Jp Jp/y 
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The 95% and 99% C.L. upper bounds, P95 and pgg, are determined by the conditions 

V{p = P95) = 0.05, V{p = P99) = 0.01 . (44) 

It is impossible to reconstruct the probabihties, PaiU'^A, Am^), from pubhshed 90% C.L. 
exclusion plots. Therefore to make estimations, we adopt the following procedure: 

1) . We assume that the distributions Pa{U^i, Am"^) have the Gaussian form. The latter is 
characterized by central value, f7^4, and by the widths Uq,. 

2) . For several Am^ the experimental groups have published the best fit values of U^^^ and 
we use them as the central values U^^ in the Gaussian distributions. 

3) . For those Am^ we find the widths of the distributions using Eqs. P^ ) and the published 
90% C.L. upper bounds on = 1/4 sin^ 26'„„. 

4) . Using definitions (g^) and Eq. (|3|), we calculate the 95% C.L. and 99% C.L. bounds, 
on the products of the matrix elements, and consequently, on Ue — mixing parameter: 
sin^ 2^eM < 4p. 

The bounds on sin^ 2^e^ are shown in Fig. § for different values of Am^ by rombs and 
triangles. Also shown, is the 90% C.L. allowed region of LSND taken from Ref. As 
follows from the figure, in the interval Am? = 0.4 — 2 eV^, some part of the LSND region 
is compatible with the Bugey and CDHS bounds at 95 — 99% C.L. . The upper limit given 
by the product of the bounds corresponds to 95% C.L. in the range of Air? < 1 eV^ and 
it approaches 99% C.L. for smaller values of Am?. 

For Am^ < 0.7 eV^, the bound from CDHS disappears and stronger restriction on 
sin^ 2^^^ follows from the atmospheric neutrino results, namely, from the up-down asymme- 
try which becomes suppressed when large admixture of exists in the heavy state 

Thus, the LSND result is consistent with bounds on mixing from other experiments at 
a few per cent level in the of the (3 + l)-scheme. 

Introduction of more than 1 sterile neutrino can enhance the LSND probability. The 
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enhancement consistent with bounds on the Ve - and - disappearance is however rather 
weak: about 20 % (see Appendix). 



4 Improving the bounds on z/^- and z/g- disappearance 

Explanation of the LSND result in the (3 + l)-scheme implies that both U^^ and ?7^4 are close 
to their present upper experimental bounds. If 1/4 consists mainly of the sterile component, 
the dominant modes of oscillations will be Ve ^ and ^ Ug. If 1/4 contains significant 
admixture of the Ur, one should expect also z/^ ^ z/^ and Ue ^ Vt oscillations [TO]. Therefore, 



further (even moderate) improvements of sensitivities of the oscillation searches with respect 
to the Bugey and CDHS sensitivities should lead to discovery of oscillations both in z/g- 
and in v^- disappearance. Negative results of these searches will exclude the oscillation 
interpretation of the LSND result in the (3 + l)-scheme. 

Let us consider possibilities to improve bounds on z/^- and z/g- oscillation disappearance 
searches in the forthcoming experiments. 



4.1 Ue- disappearance experiments 

The z/e- disappearance due to oscillations with Am^ = Am\gj^^ = (0.4 — 10) eV^ can be 
searched for in the high statistics short base-line reactor experiments of the Bugey type. 
Keeping in mind restricted power of the reactors the only possibility to significantly increase 
statistics is to increase the size of the detector placed close to the reactor. Notice that the 
absolute value of the neutrino flux is known with an accuracy about 3 - 5 %, and therefore 
to improve substantially the Bugey bound one should search for the non-averaged oscillation 
effects. Notice that for Am^ ~ 1 eV^ and E = 3 MeV the oscillation length equals Ijy ~ 7 — 8 
m which is on the border of the averaging regime. 

No experiments of this type are planned now, and it seems, the Bugey bounds on U^^ 
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will not be improved before Mini-BooNE will report its results ^0 



The reactor project can be considered if Mini-BooNE will confirm the LSND result. Such 
an experiment will allow to disentangle the (2 + 2) and (3 + l)-schemes. Indeed, in the (2 + 
2)-scheme one expects small i^e-disappearance effect which can be described by the mixing 
parameter of the order the LSND mixing parameter: sin^ 26*66 = 4(?7e3+^e4) ~ 2 sin^ ^Olsnd- 
That is, the positive effect of oscillations in the short-baseline reactor experiment will ex- 
clude the (2 + 2)-scheme. 



4.2 i/n - disappearance and KEK - front detector experiment 



Let us consider a possibility to improve bound on using KEK front detectors ^T|]. At 



KEK the high intensity z/^-fiux with energies 0.5-3 GeV (maximum at 1 GeV) is formed 
in the decay pipe of the length 200 m. Two front detectors are situated at the distance 
about 100 m from the end of the pipe. At the Fine Grained Detector (FGD) one expects to 
observe about 9200 /i - like events for 1. 10^^*^ protons on target (P.O.T.) , whereas at the 

I kton water Cherenkov detector the number of events is 30 times larger. 

One can search for the z/^-disappearance by measuring the energy distribution of the 

II - like events. Let us evaluate the sensitivity of such a study. We calculate the ratios 
of numbers of events in the energy bins i: E^^^ E^^ + with (iV*^^) and without (Nq) 
oscillations: 

Ri = ^ ■ (45) 
The ratios can be written in the following way: 

1 rEi+AE roo roc 

R^^J;^ j^^ dE,j^ dE'J{E„E'^) j^^ dE,F{E,)a{E,, E'^)P{E,) , (46) 

jj, jj, 

where f{E^,E'^) is the (muon) energy resolution function, F{Ey) is the neutrino flux at 
the front detector a^E^jEjj) is the total inclusive cross section of the u^N — > fiX 



reaction |2^. It is dominated by quasi-elastic cross-section with sub-leading contributions 
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from other reactions in the energy range of the KEK experiment (^E^"'^ < 3 GeV). In our 
calculations we use the quasi-elastic cross section only. is the normalization factor which 
equals the integral (|i6| ) with P = 1. In Eq. (^61), P{Eu) is the neutrino survival probability 
averaged over the production region: 

P{E^) = iTTTiTT / dx{L - x)-^P{E^, L - x)e — . (47) 



Np{E,) Jo 

Here 



is the effective Lorentz factor of pions which produce neutrinos with the energy E^, and 
k = 0.488, m,r and are the mass and the lifetime of the pion; L ^ 300 m is the distance 
between the beginning of the decay pipe and the front detectors, Xp is the length of the pipe. 
In Eq. (^) Np is the normalization factor which equals the same integral with P = 1. 

Since the energy resolution is rather good and the size of the energy bin is smaller than 
typical energy scale of the probability changes we omit the integration over E^j^ and E'^ and 
in the rest of integration take E'^ = E^^. 

The results of calculations of Ri are shown in the Fig. § for several values of Am? 
and mixing angles at the border of the region excluded by the CDHS experiment. The 
characteristic oscillation pattern is clearly seen for large Am? > 6 eV^. For Am^ < 1 eV^, 
we find very smooth distortion of spectrum with increasing deficit at events of low energies. 
The deficit is at most 10%. As follows from the figure, with further increase of statistics the 
KEK front detector may improve the CDHS bounds on U'^^ in the high Am^ part of the 
LSND allowed region. 



It would be interesting to estimate the possibility of water Cherenkov front detector 
where the statistics should be about 30 times higher than in FGD. The problem here is that 
the muons with i?^ > 1.6 GeV are not contained inside the detector p6| and therefore 



the higher part of muon spectrum cannot be measured. 

Let us underline that in this section we have described results of simplified estimations 
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and the detailed study of possibilities by the experimental groups are necessary. 

4.3 Mini-BooNE experiment, and the neutrino mass spectrum 

In the case of positive result of the i/^ — Vf, oscillation searches by the Mini-BooNE we 
will have strong evidence of existence of the sterile neutrino. Moreover, the Mini-BooNE 
experiment itself may allow us to disentangle the (3 + 1)- and (2 + 2)-schemes. Indeed, 
Mini-BooNE can search for the z/^ - disappearance in the range Am^ = (8 ■ 10^^ ~ 20) eV^. 
For Am^ ~ 0.6 eV^, the sensitivity to the effective mixing parameter can reach 0.15. 

In the (2 + 2)-scheme the z/^-disappearance driven by Aml^jvD has very small proba- 
bility suppressed by small admixtures of in the "wrong" pair. If is mainly in the pair 
(1/3, z/4) then 

sin^ 2^,^ = 4(f/;i + f/y . (48) 

It is at the level of mixing imphed by the LSND: sin^ 26*^^ < sin^ '^^Qlsnd- 

In contrast, in (3 + l)-scheme z/^-disappearance should be at the level of the bound from 
the CDHS experiment |]14|, as is given in Eq. (|38|) . That is, sin^ 26'^^ ~ sin^ 26'^e/?7g4 - the 
mixing parameter is enhanced by the factor \jlJl^. 

If the Mini-BooNE experiment will find Vy, — > z^e oscillation effect compatible with LSND, 
but no z/^-disappearance will be detected, the (3 + l)-scheme will be excluded, provided that 
A^i57VD = (2-4)-10^^ eV^. 

For Amj^gj^^ > 4- 10~^ eV^, the sensitivity of Mini-BooNE to z/^ disappearance is worse 
than the present bound from CDHS experiment. Much better sensitivity to z/^ disappearance 
can be reached in the BooNE experiment. In Fig. ^ we show the restrictions on the LSND 
oscillation mode from the bounds which can be achieved by the Mini-BooNE and BooNE 
experiments instead of the CDHS bound. 
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5 Phenomenology of (3 + 1) -mass scheme 



In the "minimal" version of the (3 + l)-scheme with normal mass hierarchy and small Ug — Vr 
mixing (one would expect ~ the presence of the sterile neutrino in the flavor block 
is rather small: 

E |t/siP~f/f4 + t^'4~(3-10)-10"^ (49) 

1=1,2,3 

It will be difficult to detect such admixtures in oscillations driven by Am^ and Am^^^, that 
is, in the solar, atmospheric and the long base-line experiments. Therefore, new elements 
of phenomenology of the scheme are associated mainly with the mixing of active neutrinos 
in the 4th (isolated) mass eigenstate. 



5.1 Neutrinoless Double beta decay 

If neutrinos are the Majorana particles, the double beta decay, 13(3^,^, can put bounds on 
the mixing parameters. Indeed, the (3 + l)-scheme predicts the value for the effective 
Majorana mass of the electron neutrino, m-ee, which can be accessible for the future (3l3Qy 
measurements. If the spectrum has normal mass hierarchy with the isolated state being the 
heaviest one, the dominant contribution follows from this heaviest state: 



W ^ UlJAmls^^. (50) 



The upper limit on as a function of Am|^^jy^ in the allowed region of the LSND 
result is shown in Fig. ^ The Heidelberg-Moscow present limits on mee is showed [27 . 
As follows from the figure the next generation of (3(3Qv experiments, the GENIUS (1 ton 
the EXO and CUORE [E^ will be able to test this scenario. For preferable 



version 



values Am^ ~ (1 - 2) eV^ we get from (H) ruee ~ (1 - 2) ■ lO'^ eV. 

The (3 + l)-scheme with inverted mass hierarchy in which the three active states form 
degenerate triplet with m > (0.4 — 10) eV is restricted already by present data. The bound 
on PPou rate [^] implies large mixing angle solutions of the solar neutrino problem (LMA, 
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LOW, VO) and cancellations of contributions to rriee from different mass eigenstates |30|| . 

5.2 Atmospheric neutrinos 

The survival probability for the atmospheric muon neutrinos equals: 

P,,^l- 2Ul, - AU'^,il - U% - Ul,) sin' Aml^, (51) 

where we have averaged over oscillations driven by the mass split, Amf^gj^j^. According 
to Eq. (0), the shape of the zenith angle distribution and the up-down asymmetry are 
determined by the effective mixing parameter 

sin^2e:g) = - 4 Ul,{l - U% + Ul,). (52) 

Notice that this parameter can be bigger than 1: e.g., for f/^g = 1/2 and = 2-10^^ we get 
sin^ 29atJn ~ 1.04. This can explain an appearance of unphysical values of sin^ 29atm in the 
fits of the zenith angle distributions. For small values of Amj^gj^j-,, where the CDHS bound is 
weaker or absent, can be large enough and an enhancement of sin^ can be significant. 
The presence of additional sterile state suppresses the average survival probability 

{P,,) ^ 1 - 2iU,,)' - 2^/^3(1 - Ul, - Ul,). (53) 

5.3 Nucleosynthesis bound 

Sterile neutrinos are generated in the Early Universe due to oscillations Ve ^ i^s and ^ Ug 
driven by Am\g^j~, and mixing of z/e and z/^ in the isolated 4th state. The mixing parameters 
for z/e — > and z/^ — * z/^ channels, sin 26 

es ~ 4f/g4 and sin^ 29 ~ 4f/^4, are big enough, 
so that oscillations will lead to the equilibrium concentration of Vg. Thus in the epoch of 
primordial nucleosynthesis (T ~ few MeV) one expects the presence of four neutrino species 
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in equilibrium. The current limit, Ny < 3.4, can be avoided if large enough lepton asym- 
metry (~ 10~^) existed already before oscillations driven by the Am^^jy^ became efficient 
(T ~ 10 MeV). In this case, the asymmetry suppresses oscillations in the z/g —>■ Vs and 
^ channels. Notice however, that large lepton asymmetry cannot be produced in the 
scheme itself by the oscillations to sterile neutrinos and anti-neutrinos: the mechanism ^2 



does not works. It requires the active neutrino state to be heavier than sterile state and 
very small mixing of sterile neutrinos in that heavier state. Another way to avoid the bound 
is to assume that active neutrinos were not in equilibrium before nucleosynthesis epoch 
In this case, the number of additional neutrino species can be even larger than 1. 



5.4 Supernova neutrinos 

Physics of the supernova neutrinos can be significantly modified by presence of the sterile 
component. 

The difference of the matter potentials for z/g and Vg in electrically neutral medium 
equals 



1 - 



rin 
2ne 



+ K = ^(3y;-1) + K, (54) 

where n is the total nucleon density; Ye = Ue/n is the number of electrons per nucleon; V,^ 
is the potential due to neutrino-neutrino scattering which is important in the region close 
to the neutrinosphere. 

Typical dependence of the potential Ves on distance is shown in Fig. |^ . For 



the neutrino channel the potential is negative in the central parts; it becomes zero when 
Yf. ^ 1/3, then it changes the sign, reaches maximum and then falls down with density. 
The potential changes also with time: the gradients become larger. Moreover, after several 
seconds another region with strong neutronization {Yf. < 1/3) may appear at densities below 
10^ g/cc. An appearance of this region is related to the conversion z/g — > z/^ in the outer 
resonance (see later). As a consequence, the reaction Ue + p + n will dominate over 
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Ue + n ^ e +p, thus producing the neutron reach region. It was suggested that in this 
region the nucleosynthesis of heavy elements may take place. 

For the z/^ — channel the difference of potentials equals 

V,s = ^{l-Y^) + Vl. (55) 

Crossings of the Ug- and Ug- levels are determined by the condition: 

^%§^-V„ = 0. (56) 

For the anti-neutrino channels Ves or the first term should be taken with opposite sign. 
Similar equation should be written for the potential V^^. 

In Fig. ^the points of the resonances are shown as crossings of the Ves with the horizontal 
lines given by ~ /S.m? /2E (lines below V = Q correspond to the anti- neutrino channel). 

From (0) and (^) we find the resonance densities (for = 0): 

_ Am|gjv^cos2g V2 

2E Gp{3Ye-l)- 

Using the potentials (|5^ ) , (|20|) and resonance conditions ( ^61) we can construct the level 
crossing scheme. The scheme for the (3 + 1) spectrum with normal mass hierarchy is shown 
in fig. |[ We use the flavor basis {ug, i/g, z/*, u*), where z/*, u* are the states which diagonalize 
the Ufj, - Vr sub-matrix of the total mass matrix. This rotation allows us to remove large 
mixing and it does not change the physics since the produced and Vr fluxes are practically 
identical and they can not be distinguished at the detection. 

The adiabaticity condition in the z/g — z/^ resonances can be written as 7 >> 1, where 
the adiabaticity parameter equals 

W f.. r 



^ = 47^2^14^ sin^ 2^ 



L\ 3 



(5^ 



Here ly = Ye/{dYe/dx), In = n/ (dn/dx), = AttE/ Am^ , Iq = \/2mM / Gpp are the vacuum 
oscillation length and the refraction length. 
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In the case of normal mass hierarchy in the neutrino channels there are two crossings 
of the Vs and z/g levels determined by the condition (|^) (see Fig. One (outer) crossing 
occurs at relatively low density pr ~ lQ'^g/crn?{^m\gj^j^/leV'^) [|^, when Ye substantially 
deviates from 1/3 {Y(. ~ 0.5) (see Fig. |]). It is characterized by the mixing parameter 
sin^ 26'es ~ 4?7g4 (4f/g4Cos/5^ in general). This resonance is highly adiabatic: the resonance 
density is relatively small, so that the refraction length is large. Also the gradients of the 
density and electron number are rather small. For E < MeV, and sin^ 26 = A - 10~^ we 
find that 7 > 0(10) in the relevant range /S.m\g^j^ = (0.4 -h 10) eV^. 

Thus, in the outer resonance the Uf, converts almost completely to Us with survival 
probability 

Pee ~ U^^ < (2 - 3) ■ 10-2 . (59) 

If Amj^gj^^ > 2 eV^, the Ue-fiux is absent above the layers with p ~ (10^ — 10^)g/cc, that 
is, in the region where r-processes can lead to synthesis of heavy elements |^ . An absence 
of the z/e-flux enhances an efficiency of the r-processes. For Am^^jv^ < 1 eV^ the effect of 
conversion on the r-processes is weak. 

The second (inner) crossing occurs at much higher densities p ~ 10^^ g/cm^ in the 
layers with significant neutronization when 

Ye ^ 1/3 

(see Fig. Here the neutrino- neutrino scattering gives substantial contribution to Ves 
which shifts the resonance. The level crossing occurs mainly due to change of the chemical 
composition (the ratio n„/ne = rin/np). Taking n ■ rriN = 10^° g/cc, (mjv is the nucleon 
mass), ly = 3 km, E = 10 MeV, Am^ = 10 eV^ and sin^ 26^65 = 0.1 we get for the 
adiabaticity parameter 7 ~ 0.3, that is, the adiabaticity is broken. For Am^ > 30 eV^ 
the adiabaticity could be satisfied. This leads to disappearance of the z/g-fiux. Therefore 
observation of the anti-neutrino signal from SN1987A excludes such a possibility thus putting 
the upper bound on Am^. (Notice that the Ue-Hux regenerates in the outer resonance 
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destroying conditions for the r-processes). At later time, gradients of density and Ye become 
larger, in particular, ly ~ 0.3 km, so that the adiabaticity breaks down even for large Am^. 

For Am^ < 2 eV^ the adiabaticity is strongly broken for neutrinos of the detectable 
energies, E > 5 MeV, during all the time of the burst and the effect of the inner resonance 
can be neglected. 

Let us consider properties of the neutrino bursts arriving at the Earth detectors. 

1) . The generic prediction of the scheme is disappearance of the "neutronization peak" 
due to z/g — Ug conversion in the outer resonance (pQ]) . The peak should be seen neither 
in z/g nor in z/^/z/,-. Observation of the peak in the z/g luminosity at the initial stage will 
exclude the scheme with normal mass hierarchy. Only if the inner resonance is adiabatic 
the scheme can survive. 

2) . During the cooling stage one expects hard z/g-spectrum which coincides with original 
spectrum of and/or z/,-: -F(z/g) oc F^{u^). The absolute value of the z/g-flux depends on 
the mixing in the flavour block. If U^^ > 10~^, so that the level crossing related to Am^^,^ is 
adiabatic, one gets F{ve) = F^{i^ti) independently on solution of the solar neutrino problem. 
No Earth matter effect should be observed. If f/gg < 10~^, the level crossing is strongly 
non-adiabatic and the z/g-signal depends on properties of the low density resonance related 
to the AttIq. For the LMA solution one gets: F{ve) = cos^ 9qF^{v^). The flux can be larger 
in the case of SMA (small mixing angle) solution depending on the level of adiabaticity. 
Intermediate situation is possible, if t/g3 is in the interval 10^^ — 10^^. 

In the antineutrino channels, the z/^-level does not cross Ve level but it has two crossings 
with z/* and v* levels which can be both adiabatic. Still the presence of the sterile neutrino 
will modify the z/g signal. The z/g propagation is adiabatic so that z/g — > z/i. If SMA solution 
is correct, one gets vi ~ z/g, that is z/g-flux is practically unchanged. For LMA (or LOW) 
the soft part of the spectrum will be suppressed: 

F(z^g) = cos2^0F°(z^g) (60) 
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and the total Ve flux will depend on characteristics of Vs — v^/vr resonances. At least 
one of these resonances should be adiabatic since the mixing of Vs in the active block is 
I]f=i V^i > f^e4 + Uf^i and the later is large enough. If crossing with i/* is adiabatic, 
then z/* — > Vg and the Vf. flux will be as in Eq. (|60D. If the resonance is completely 
non-adiabatic, then i/* 1^2, and since 1^2 has a component sin^ 6'q of the u^-^ux we get 
F(z/e) = cos^ 6'0F°(z/e) + siu^ 6*© F° ( z/^ ) . In these cases also the Earth matter effect should 
be observed. 

Some combination of the z/^ and Ur fluxes, will be converted to the sterile neutrino flux. 

In the case of inverted mass hierarchy, when the isolated fourth state is the lightest one, 
the crossings of Ue and Ug levels appear in the antineutrino plane. As the result, originally 
produced z/g will be converted to Ug- For Arnj^gj^^ > 2 eV^ this will dump the efficiency of 
the r-processes. 

The z>e flux at the Earth will be formed in the conversion of z/^ and Ur to z/g due to the 
mixing in the flavor block: more precisely, in the resonance associated with AttIq. In the 
case of SMA solution this conversion is practically absent, so that no z/g-flux is expected 
at the Earth in contradiction with observations of the z/g-burst from SN1987A. Therefore 
the SN1987A data require one of the large mixing solutions. In the case of LMA or LOW 
solutions one expects -F(z/e) ~ sin^ 6'0F°(t'^). Properties of the z/g-flux are determined by 
mixing in the flavour block ||30| . 

6 Discussions and conclusions 

1. Simultaneous explanation of the LSND result as well as the solar and atmospheric neu- 
trino data in terms of neutrino oscillations requires introduction of additional (sterile) neu- 
trino. In the favored (2 + 2) - scheme the sterile component should be involved signiflcantly 
in the conversion of the solar and atmospheric neutrinos. 

2. Recent experimental results give an indication that both in solar and in atmospheric 
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neutrinos the oscillation channels to active component dominate. We show that in the (2 
+ 2)-scheme the equahty rjf^'^ + rjf'^ — 1 should be fulfilled. The statistically significant 
deviation from this equality, in particular, a proof that rjs = rif^^ + rif"^ < 1 will lead to 
rejection of the (2 + 2)-scheme. 

In the case of (3 + l)-scheme rff^ can take any value between ^ and 1 + sin^ Oatm- 
The generic prediction of the scheme is that 77!"" < rif^. 

Determination of r]s from analysis of the solar and atmospheric neutrino data will allow 
us to disentangle the (2 + 2) and (3 + l)-schemes. 

We show that the LSND result can be reconciled with the bounds on z/g— (Bugey) and 
v^— (CDHS) disappearance at 95 % - 99% C.L.. The best agreement can be achieved at 
Am2 = (l-2) eV^. 

4. The key consequences of the (3 + l)-schcmc which explains the LSND result are that 
the probabihties of the ^ and ^ oscillations with Am^ = Am^^^^, should be 
at the level of present upper bounds, so that moderate improvements of the bounds should 
lead to discovery of oscillations in these channels. We have estimated the possibihty of the 
KEK - front detectors experiments to improve the bound on i/^ disappearance. We show 
that for Am^ > 6 eV^ searches of the distortion of the energy spectrum of the /x-like events 
may have better sensitivity to oscillations than CDHS has. The Mini-BooNE experiment 
can improve the CDHS bound for Am^ < 0.6 eV^ and the BooNE experiment will have 
better sensitivity than CDHS for Am^ < 1 eV^. 

At the moment, no experiment is planned to improve the Bugey bound on i/e-disappearance 
in the range Am^ = ■^^isiVD- The relevant experiment could be the one with large size 
detector situated close to the reactor. Such an experiment will deserve serious consideration 
if Mini-BooNE confirms the LSND result or/and some indications of the sterile neutrino 
involvement are seen in the solar and atmospheric neutrino experiments. 

5. In the (3 + l)-scheme the effective Majorana mass of the electron neutrino can be 
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about 2 ■ 10 ^ eV for the normal mass hierarchy and at the level of the upper experimental 
bound in the case of inverted mass hierarchy. 

6. The scheme can reproduce an "unphysical" value of the effective mixing parameter: 
sin^ 29'^^^ > 1 in the two neutrino analysis of the atmospheric neutrino data. 

7. The scheme predicts that a sterile neutrino was in the thermal equilibrium in the 
epoch of primordial nucleosynthesis unless "primordial" lepton asymmetry suppressed the 
oscillations of active to sterile neutrinos. 

8. Presence of the sterile neutrino in the (3 + l)-scheme can substantially influence 
the physics of supernova neutrinos. The (3 + l)-scheme with normal mass hierarchy and 
Am^ > 2 eV^ can realize mechanism of disappearance of the Ug-finx from the region with 
densities (10^ — 10®) g/cc, thus creating necessary conditions for the nucleosynthesis of the 
heavy elements in the r-processes (see however |}46[ ) . 

The scheme predicts disappearance of the neutronization peak and hard u^- spectrum 
during the cooling stage. The presence of sterile neutrino can manifest itself also as sup- 
pression of the Ue flux without change of the spectrum. 

Let us outline further developments. 

1). Before Mini-BooNE result, the progress will be related to further searches for the 
sterile neutrino components in the solar and atmospheric neutrinos. The experimental 
groups are requested to get bounds on rj'^'^"' and i]"*'" parameters from the one Am^ fit of 
their results. 

Important bounds will come from the SNO experiment: namely, from measurements of 
the CC event rate and comparison of this rate with the neutrino-electron scattering rate, 
from searches for the Earth regeneration effects and the distortion of the spectrum, from 
measurements of the NC-event rates and the double ratio [A^C]/[CC]. 

As far as the atmospheric neutrino data are concerned, further accumulation of data on 
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the zenith angle distributions and neutral current interactions will strengthen the bound on 
^atm^ Here the distributions of the upward going muons and the events in the NC - enriched 
samples, as well as measurements of 7r° events rate are of special interest. 

Still it is difficult to expect that the bound rjs < 1 will be established at the high 
confidence level. 

2). The Mini-BooNE experiment will give the key result which will determine further 
developments. In the case of negative result, still searches for the sterile neutrinos will be 
continued with some other motivation. The positive result will give strong confirmation of 
existence of the sterile neutrino, and identification of correct four (or more) states scheme 
will be the major issue. 

Notice that values of the sin^ 29e^ at the upper side of the LSND allowed region will 
favour the (2 + 2)-scheme. Also restrictions on Am^ (say excluding of Am^ > 1 eV^) will 
have important implications. 

Mini-BooNE itself and BooNE will search for also z/^-oscillation disappearance. Detec- 
tion of the disappearance will favor the (3 + l)-scheme. 



The ORLAND experiment |^9[ will have sensitivity for z/g oscillations which will 

cover the LSND region. 

3). Later progress will be related to searches for the Ur— appearance and z/^— disap- 
pearance in the long baseline experiments. OPERA ||5^ will allow one to improve bound 
on the sterile neutrinos. MINOS experiment |^ will provide a decisive check of presence 
of sterile neutrinos by measurements of NC/CC ratio. 

Special high statistics reactor experiment to search for the z/g-disappearance can be 
considered. 

Future neutrinoless double beta decays searches can give an additional check of the 
model being especially sensitive to scheme with inverted mass hierarchy. 
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Appendix. LSND and additional sterile neutrinos 



Let us make several remarks on possibility to get larger predicted LSND effect by introducing 
additional sterile neutrino v'^. In general an introduction of additional neutrino state can 
enhance the LSND probability and/or influence (relax) bounds from the laboratory and 
atmospheric neutrino experiments on the relevant mixing parameters Uei and 

1) . To influence the prediction for LSND oscillation probability, the 5th mass eigenstate 
(being mainly sterile neutrino) should have non-zero mixing with both v^, and v^. If z/5 
has admixture of z/g (or z/^) only, it will open new channels for Ve-{i-'fi-) disappearance, 
thus enhancing the bound on U^^ (or f/^4) from Bugey (CDHS) experiment. At the same 
time, it will be no additional contribution to the LSND probability in comparation with the 
(3 + l)-scheme. As the result, one gets stronger bound on possible value of sin^26'e/i, 

sin^ 20,, < 4 f/L,e,(A^')^CDH5(Am2)- (61) 

2) . The sensitivities of the Bugey, LSND and CDHS results to Am^ have the following 
hierarchy: /S.m\^g^y < Am^^^^, < /S.m'^j^^g. Then depending on value of m\ one can get 
different situations: 

(i) m\ < Am^^ggj^: neither bounds on nor the LSND probability are changed. 

(ii) /S.m\^g^y < m\ < Am^^^^,: No additional contribution to the LSND appears. At 
the same time, new channel will be open for the z/g-disappearance which leads to stronger 
bound on U^j^. As a consequence, the LSND probability will be suppressed. 

(iii) m\ > /S.m\g^j^: for mixing angles relevant for the LSND at this condition both 
bounds on will be modified and the LSND will get an additional contribution from the 
second Am^. 

The mass should not be very close to 777,4. The splitting should be resolved by the 
LSND experiment, that is: Am54 >> 10~^ eV^. Following the analysis of Ref. [Q, we find 
that in the degenerate case, Am54 ^ 10~^ eV^, the depths of the z/^- and z/g-oscillations 
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equal 

sm^2e^^ = 4r]^{l-7]f,), (62) 

Sin2 2^ee = 477e(l-r/e), (63) 

where rja = J2i=4,5 U^i: a = e, fi. The depth of the z/^ — > Ue channel driven by Aml^jvD is 
sin^ 2eLSND = 4 E U^^Ul < 4 ^ [/J ^ [/^ = r^,r^^ = i sin^ 29,, x sin^ 2^^^ , (64) 

1=4,5 i=4,5 i=4,5 



where we have used the Schwartz inequality. Comparing result in (0) with the one in ([39| ) 
we see that the LSND probability in the scheme with additional neutrino is smaller than in 
the (3 + l)-scheme. 

As an example, we have considered the neutrino mass spectrum with (i) Am|. ^ (1 — 
2) eV^, where index i = 1,2,3 enumerates the eigenstates from the flavour block, (ii) 
Amgj > 8 eV^ and (iii) non-zero admixtures, 11,5 and of the and i/^ in the 5th state. 
So, both necessary conditions discussed above are satisfied. Oscillations associated with the 
5th state are averaged in the Bugey and LSND experiments and they are at least partly 
averaged in the CDHS experiment. 

In the mass range, Am|j > 8 eV^, the CDHS bound is rather weak and stronger bound 
on the follows from the CCFR experiment This CCFR bound on is still a 

factor of two weaker than the CDHS bound on at Am^ = (1 — 2) eV^. Taking also 
f/fg ~ f/g4 ~ U]j^g,y/2 we have found that the probability of z/^ v, oscillations can be 
enhanced by a factor 1.25 at most with respect to result of the (3 + l)-scheme (see Eq. (0)). 

Notice also that introduction of the second sterile neutrino can further aggravate a 
situation with primordial nucleosynthesis. 
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Note added 



After this work was accomplished the paper [52] has appeared in which the possibihty to 
constraint the t]^**" parameter to be down to 0.3 using the energy distribution of the CC 
and NC events at MINOS was discussed. 
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Figure 1: The neutrino mass and flavor spectrum in the (3 + 1) scheme with normal mass 
hierarchy. The boxes correspond to the mass eigenstates. The sizes of different regions in 
the boxes determine flavors of the eigenstates: \Uai\'^ , a = e, ii,t, s, i = 1, 2, 3, 4. 
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Figure 2: The LSND allowed region (shadowed) and the bounds on the oscillation pa- 
rameters in the (3 + 1) - scheme from the z/g— and z/^— disappearance experiments. The 
lines show the limits obtained as the products of the 90 % CL. upper bounds on |?7e4p and 
|f/;_j4p. Solid line is product of the bounds from Bugey and CDHS (above the arrow) and 
is product of the bounds from Bugey and atmospheric neutrinos (below the arrow). Also 
shown are the bounds obtained as the product of bound from Bugey and expected bound 
from Mini-BooNE (dashed line) and BooNE (dotted line) experiments. The triangles and 
rombs show respectively the 95 and 99 % CL. bounds obtained in this paper (see Eq. ^31) . 
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Figure 3: The upper bound on the effective Majorana mass of the electron neutrino rriee 
as the function of ^rn^sND (3 + 1) scheme with normal mass hierarchy. The bound 

corresponds to the upper limit on |C/e4p from the Bugey experiment. Also shown are the 
present and future limits on niee from the /3/3oi/ - experiments. 
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Spectre II distortion of KEK Front Detector 
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Figure 4: Distortion of the energy spectrum of the //-hke events due to oscillations in the 
KEK - front detector experiment. Shown are the ratios of the predicted numbers of events 
with and without oscillations for different values of the oscillation parameters. 
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Figure 5: The matter potentials for the Ue — Vg system in the supernova as the function 
of distance from the center of the star. Sohd hne shows the potential without neutrino 
background and the dashed line - with the background. The horizontal lines correspond to 
±Amj^gj^jj/2Ei, for E^, = 10 MeV and Am^ — 10 eV^; the positive (negative) sign refers to 
neutrino (antineutrino) channel. Crossings of the potentials and the lines give the resonance 
points. 
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Figure 6: The level crossing diagram for the (3 + 1) scheme with normal mass hierarchy in 
supernova. Solid lines show the dependence of the eigenstates of the effective Hamiltonian 
on the matter density. Dashed lines show the dependence of the flavor states. The semi- 
plane with n > (n < 0) corresponds to the neutrino (antineutrino) channels. Only the 
outer resonance v^. — Vg is shown. 



